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Stability of Magnetically Trapped Bose-Einstein Condensates 
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According to the adiabatic approximation atoms moving in a magnetic trap keep their mag- 
netic states. We investigate the validity of this approximation for quantum condensates, where 
a change of field's direction generates effective interactions between hyperfine angular momen- 
tum states. Condensates in general traps are found to be stable because they are confined in 
the vicinity of the trap center. A decay of a condensate is observable in a trap with extremely 
large field gradient. 
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1. Introduction 

Bose-Einstein condensate (BEC) of ultracold atoms 
has been studied furiously since its experimental achieve- 
ments. 1-3 ' It still provides us a lot of interesting features. 
Recently a quite small condensate on a microelectronic 
chip was realized, where modest electric currents pro- 
duced large magnetic field gradients and a microscopic 
trap was formed on a chip. 4 ) A quantized vortex in a 
condensate is also one of the interesting phenomena. A 
method of creating a vortex topologically through in- 
verting the bias field along the trap axis was suggested 5 ' 
and experimentally confirmed. 6 ' Since the field gradi- 
ents used in these experiments are large as compared to 
the bias fields, spin of a condensed atom may flip in the 
trap. 

According to the adiabatic approximation, an ultra- 
cold atom moving in a spatially changing magnetic field 
keeps the relative direction of its magnetic moment to the 
field. Atoms are classified into three states by the sign of 
Zeeman energy: weak- field-seeking state (wfss), neutral 
state (ns) and strong-field-seeking state (sfss). In a mag- 
netic field which has a minimum point of the magnitude 
in space, only the atoms in wfss can be trapped. 

The adiabatic approximation requires that the field's 
direction 6 which trapped atom feels changes slower than 
the precession of the magnetic moment: 



dO 

~~ 77 ^ ^Larmor- 

at 
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This adiabatic condition is violated in regions of very 
small magnetic fields. The motion of spin in a trap with 
no bias field has been calculated classically. 7 ' It was 
found that the size of such non-adiabatic region is very 
small and the trapped atoms satisfies the adiabatic con- 
dition. 

In this paper we consider quantum condensates. When 
we describe the Gross-Pitaevskii (GP) equation with or- 
der parameters corresponding to wfss, ns and s/ss, 8 ' cou- 
pling terms between these states appear in the equation, 
besides nonlinear terms. The adiabatic approximation 
requires these coupling terms to be negligible. We evalu- 
ate the effects of these coupling and investigate the limit 
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of the application of the adiabatic approximation by cal- 
culating the lifetime of a condensate. 

We treat alkali atoms with a hyperfine state |F| = 1 
such as 23 Na and 87 Rb. In the next section we briefly 
introduce the treatment of BEC with internal degrees of 
freedom, we represent the GP equation for wfss, ns and 
sfss, and derive the expression of the decay rate. In Sec. 3 
the lifetime of a condensate is calculated numerically. We 
discuss about its dependence on field geometries. The 
effect of singularity of the condensate and population of 
the trapped atoms are also referred to. The final section 
is devoted to conclusions. 

2. Bose-Einstein Condensation in a Magnetic 
Field 

2.1 Gross-Pitaevskii equation o/|F| = 1 BEC 

Let F be the angular momentum operator with hy- 
perfine state |F| = 1. The eigenvalues of F z are 1,0 and 
— 1. The corresponding eigenstates are |1), |0) and | — 1), 
which satisfy F z \m) = m\m). In this basis, the order 
parameter |\&) is expanded as 



(2) 



m=0,±l 



It is convenient in a certain case to take another basis set 
of \x), \y) and \z), which satisfy F a \a) =0 (a = x, y, z). 
|W) is represented in this basis as 



(3) 



Clearly, the state \z) corresponds to |0) and other states 
satisfy the relations, given by 



\±l)= T ^{\x)±i\v)). 



(4) 



The Hamiltonian, which should be invariant under 
spin space rotation and gauge transformation except for 
Zeeman term in the xyz-representation, is described as 



H 
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where 7 M ~ — ^b/2 is the gyromagnetic ratio and tp^, tp a 
are the field operators for a Bose particle in the \a) state. 
The coupling constants are given by 



9i 
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The scattering lengths for 87 Rb are clq — 5.82 nm and 
&2 = 5.66 nm. 9 ' From the equation of motion for ip a , 
one obtains the Gross-Pitaevskii equation: 



dt 



(7) 



where the ground state average (tp a ) is replaced by \& a 
and approximations such as (ip^ipip) <~ VJ/* \Er vj> are m _ 
troduced. The GP equation for ^+1^,-1 is obtained by 
transforming ^ x , y ,z to W+i^-i in cq. (7). When the 
field is homogeneous in space and parallel to z-axis, the 
Zeeman term yields no coupling in the representation of 
^+1,0, -1 and is simply described as \'j ll \Bm^> m . 

2.2 Magnetically trapped BEC 

Let us consider a condensate in a magnetic trap. The 
magnetic field must have a minimum point of its mag- 
nitude to trap atoms. In most instances the magnetic 
traps used in laboratories, for example a Ioffe-Prichard 
trap, make spatial change of the field's direction as well 
as the magnitude. This change causes spin- flips. 

We assume a Ioffe-Prichard trap with a quadrapolc 
magnetic field 



B±(r) = 



^B±(r) cos(-0) N 
B±(r) sin(— 4>) 




(8) 



where r = (r, </>, z) is the cylindrical coordinates. The 
magnitude of the quadrapolc field is proportional to r 
near the trap center: B±(r) ~ B'r at r ~ 0. Assuming 
that the trap size along to z-axis is much larger than its 
perpendicular size, we approximate the z-component of 
the magnetic field to be uniform: B = B z. This bias 
field B should be strong enough to satisfy the condition 
1 7ju \Bq ^S> ksT so that the thermal fluctuation does not 
dominate in spin-flip processes. 

In such a spatially changing magnetic field, we chose 
the direction of the local magnetic field as a quantum 
axis 



I = (sin [3 cos 0, — sin (3 sin <fi, cos/3). 



(9) 



Here (3 is an angle between z-axis and the magnetic field, 
given by 



(3 — arctan ( — 77— - I : arctan 
V B 



(10) 



where R = Bq/B' is a characteristic length of the field 
gradient. The state |tu(r)) corresponding to wfss, in 
which a hypcrfinc spin keeps its direction parallel to the 
magnetic field, satisfies 



(F-0Kr)) = Kr)). 



(11) 



Similarly, the states |n(r)) and |s(r)) representing ns and 
sfss satisfy 

(F-0|n(r))=0, (12) 
(F.f)Kr)) = -Kr)>. (13) 
These states \i(r))(i = w, n, s) are written with \a) as 
|/(r)) = >; U'Ja). (14) 
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By the analogy with the relation between |1, 0, — 1) and 
\x,y,z), the transform matrix W is defined with the 
field's direction I (and therefore depends on space) as 



W = e" 



V 72^-^) 



(15) 



where m,n are real vectors satisfying m x n = I and 
(/, m, h) forms a triad. The transformation has an ar- 
bitrary phase, and we choose it as —<j> so that the state 
|u>(r)) has no singularity at the trap center ((3 = 0). 
Then the winding numbers of the states |n(r)) and |s(r)) 
are —1 and —2, respectively. 

The order parameter |*) can be expanded with the 
states \i(r)) as 



l*> = 



-w.n.s 



*i(r)l*(r)>, 



(16) 



The transformation between ^ a and is given by 

^ a = {U' ic f^ i =U ai ^ i , (17) 

where the subscripts a,i denote x,y,z and w,n,s, re- 
spectively, and repeated index is summed. Hereafter we 
use these characters in this notation. 

Substituting the transformation eq.(17) into eq. (7), 
we can construct the GP equation for as 



.91 (*;*j)*< 

\ lf t\B(r)(ie af}x huLUp^ 3 . 



(18) 



The first term in the right hand side describes the ki- 
netic energy, which yields the additional terms W(\7U) ■ 
(V*) +W(W 2 U)^. These extra terms come from the 
spatial change of the field's direction and cause coupling. 
Since we have chosen the quantum axis in the direction 
of the magnetic field, the Zeeman term does not yields 
coupling term and the matrix elements are simplified as 



—\e a p\i\U} a Upj = Sij(5 W i - 6 si). 
The matrix U T U also have a simple form 




(19) 



(20) 
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The symmetry of the field suggests that \E', is inde- 
pendent of z. Assuming a cylindrical symmetry of the 
condensate, we divide r and ^-dependence of the order 
parameter as 



* 4 (r) = .f,(r,t)e i; 



(21) 



From eq.(18) ^-dependence of each order parameter 
should satisfy the relation m w = m n = m s = m. This 
means that wfss with winding number m can couple only 
with ns with singularity m — 1 and sfss with singularity 
TO — 2. 

When the distribution size of the trapped atoms is 
much smaller than R, the Zeeman term is approximated 
as a harmonic potential 



l7p|SW = l7pl V / ^M 2 +5 2 



(22) 



It is convenient to scale lengths and energies by the har- 
monic oscillator length oho and the energy fiwno, where 



oho = 



who 




(23) 



(24) 



The length a H o is m the order of the distribution size 
of the condensate. When we use r and / as dimension- 
less variables (scaling r/ano — > r,/4o — * /)j the GP 
equation is rewritten in a simple form, given by 



i dfi 



1 



v ijfj ~ frfi+9i(fjfj)fi 



who dt 2 

+ B(r)(d wi - S si )f h 

where B = /hwno, A = m/^ho, 9i = W a HO^no- 
The matrix T> comes from the kinetic energy term 
V 2 (U^f) and consists of the functions of r,<f> and the 
differential operators with r. 

2.3 Lifetime of weak- field- seeking state 

To calculate the lifetime of trapped condensate, we 
consider the situation that the only atoms in wfss exist 
in the trap at the time t = and populations of ns and 
sfss increase with time. 

Assuming \f w \ » \f n \, |/ s |, wfss is regarded as a sta- 
tionary state. When we trap N atoms par oho along 
z-axis and normalize f w (r) by 



/>oo 

2vr / ralr\f w (r)\ 2 = 1, 
Jo 



wfss obeys the equation 



frfu 



-^D ww f w + ~ 9l N fJ + B(r)f w , 



(26) 



(27) 



where V ww is given by 



=©o(m) + 



2m(l - cos/3) 



f3' 2 -7 + 8cos/3-cos2/i 
~1T + 4T 2 ' 

cP 1 to 2 
dr 2 r dr r 2 

Similarly, the equation for ns is also simplified as 

i 



(28) 



(29) 



It) = [ P °( m - + ' ( 3 °) 

who Ci/ ^ 



(cos /3 - 2) sin /3 + r/3' + r 2 /3" 

Vl^ 2 



(31) 



Since the kinetic energy term T> nn has no singularity at 
r = except for T> (m — 1), we neglect the additional 
terms in V nn . We also neglect spin- flips to sfss because 
wfss and sfss do not interact at r = 0. 

It is very useful to expand /„ with the (to — l)th Bessel 
function J m _i(fcr) as 

/oo />oo 
dwe^ 4 / kdkJ m ^(kr) Jn (k,Lu). (32) 
-oo JO 

Substituting this form into eq. (30) and using the Bessel 
function's property: 

pOO POO 

g(x) = / ydyj m (yx) / zdzJ m (yz)g(z), (33) 
Jo Jo 

we obtain a description of /„ (r, t) , which should satisfy 
the initial condition f n (r, 0) = 0, as 



/n(r,t) 



Jo 



kdkJ m - 1 (kr)F w (k) 



exp[i(/x - fc 2 /2)w H o*] - 1 



fi-k 2 /2 



(34) 



where 



F w (k) = l [ rdrJ m ^(kr)V nw f w (r). (35) 
^ Jo 

Then the lifetime of the condensate, or the rate of atoms 
which flip to ns par unit time, is given by 



i = 2tt r 



rdr\f n (r,t)f 



(36) 



47T WHO 



^(V 2 A) 



where the concerning time scale is much longer than 
the period of Larmor precession and the approximation 
(sinwi/w) 2 <~ 7rf<5(w) is applied. 

The function F w (k) depends on the distributions of 
wfss and ns, and the lifetime r does too. Since the atoms 
in ns do not feel Zeeman potential, they are in an eigen- 
momentum state J m -i(kr) with energy fi, which corre- 
sponds to the chemical potential of wfss. We show in 
the next section that this energy ft strongly affects the 



4 



Yuki Kawaguchi and Tctsuo Ohmi 



coupling. 

3. Calculation and Discussion 

Here we introduce a new parameter 6 = R/ano, which 
indicates the changing rate of the field in the region of 
w/ss's distribution. From the definition of R and ano: b 2 
corresponds to the ratio of the Zeeman energy at the trap 
center to the harmonic oscillator energy. We investigate 
two cases of b » 1 and b ~ 1, respectively. In the most 
cases of experiments b is much larger than 1, for instance 
b = 51 for 87 Rb in a trap of B = 1 G, B' = 300 G/cm 
and b = 10 in B = 1 G, B' = 8 kG/cm. 

3.1 Distribution in the small field gradient b ^> 1 

When the distribution of the condensate is restricted 
within a narrow region at the trap center, b 3> 1, the ap- 
proximation eq. (22) is valid. In this limit the magnetic 
field is almost parallel to z-axis ({3 <C 1). The operator 
V ww and V nw can also be expanded with 1/6. Neglecting 
the higher order of 1/b, these terms are simplified as 



V ww = T> (m), 



V — 



_ V2 ( d 



b \dr 



m 
r 



In this case, F w (k) is described with a simple form: 
k 



F w {k) 



V2b 



I 

Jo 



rdrJ m (kr)f w (r). 



(37) 
(38) 

(39) 



It should be noted that using the symmetry of the Bessel 
function J_ m (fcr) = (— l)" 1 J m (kr), the lifetime does not 
depend on a sign of the winding number m. The life- 
times of wfss with the winding number 1 and —1 arc 
same although one couples with ns with no singularity 
and the other couples with ns with a double quantized 
vortex. The stability of wfss is not affected by the wind- 
ing number of ns since once an atom flips to ns the atom 
is no more trapped. 

When the condensate density is low, the non-linear 
term giN f^ in eq.(27) is negligible. The analytical so- 
lution is given by 



Mr) = 



r \m\ e -r 2 /2 



p, = b 2 + \m\ + 1 ~ b 2 . 



Then F w (k) and r 1 are as follows: 



F w (k) 



j 1 1. k \m\ e -kV2 

y/27r\m\\ b 

(2/i)l™l _ 2a 
to ! 



(40) 
(41) 

(42) 
(43) 



where the approximation fi ~ b 2 is used. The lifetime is 
rapidly increasing function of ft and therefore of b roughly 
as exp(26 2 ). When the trap geometry is b = 10, the con- 
densate has such a long lifetime as r ~ e 200 . Although 
the singularity of the condensate shortens the lifetime, its 
effect is much smaller than the exponential factor. The 
condensate with a few winding number is stable after all. 



When the population of trapped atoms is large, it is 
necessary to solve eq. (27) numerically. The distribu- 
tion of wfss expand outward as No increases. The 6- 
dependence of the lifetime is roughly same as that of 
analytical calculation: the lifetime rapidly increases as 
b increases. To investigate its dependence upon N , we 
have particularly calculated in the case of b = 7. Sup- 
pose B — 1 G, this condition implies B' — 1.6 T/cm and 
who/2tt = 14 kHz for 87 Rb. Such a large field gradient 
may cause other experimental difficulty in achievement 
of BEC in the trap. 

Figure 1 shows our result for 1/r as a function of 
No/auo, which is a number of the trapped atoms par 
unit length. We have investigated condensates with 
m = 0, 1,2. The lifetime of wfss with to = is much 
longer than others and 1/r ~ in Fig. 1. The lifetime 
becomes shorter as N increases in the graph. The de- 
pendence of the lifetime upon No is not monotonic in 
regions of larger No (out of the figure). The interaction 
between wfss and ns is affected by the distribution of 
f w (r), and the condensates never decay for certain dis- 
tributions. An example for such point is 3300 atoms//zm 
in double quantized vortex state. 

Let us compare this lifetime with the rate of three- 
body recombination. In this process two spin-polarized 
atoms form a molecule and free from a trap. The decay 
rate is given by 



dn , 

Tt = - Ln > 



(44) 



where n is a condensate density and the recombina- 
tion coefficient for 87 Rb is L = 4 x 10~ 30 cm 6 /s. 10) 
The lifetime defined by eq. (44) is r roc <~ 1 sec at 
n ~ 5 x 10 14 cm~ 3 . Such density corresponds to 
No/auo ~ 50 ^irn -1 in Fig. 1, where r is much longer 
than 1 sec at this point. 

From these considerations, as far as b 3> 1 the conden- 
sate is found to be stable and the adiabatic approxima- 
tion is valid. It is also mentioned that the condensates 
achieved in laboratories were considerably stable with 
this decay mechanism. 

3. 2 Distribution over the large field gradient b ~ 1 

Next we consider the case of b ~ 1, where atoms mov- 
ing in a condensate feel large change of the magnetic 
field and the condensate decays faster. In this case it 
is impossible to use eqs. (22), (39). The 6-depcndcncc 
of the lifetime is exp(6 2 ) approximately, and comparable 
to T roc at b <~ 3. This is the limit that the adiabatic 
approximation is valid. 

Figure 2 shows the result in a trap of 6 = 2, B = 0.1 G 
(i.e. B' — 6.1 kG/cm) as functions of iVo/ano- The 
distribution size is similar as that of Fig. 1, and r roc <~ 
1 sec at No/auo ~ 50 /xm _1 . The lifetime at this point 
is much smaller than 1 sec. The decay of condensates 
is observable in this trap. It should be noted in Fig. 2 
that against our intuition the lifetime becomes longer 
as iVo increases. The reason is that trapping a large 
number of atoms brings an increase in the energy ft as 
well as a change of distribution f w (r) . The fc-dependence 
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m = 1 



m = 2 
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50 100 150 200 250 300 350 400 450 500 
N /a HO [nm- 1 ] 

Fig. 1. The lifetime in trap of b = 7; The field is B = 1 G, 
B' = 1.58 T/cm and ajjo = 90 nm. The condensate of m = 
is much stable and 1/r = in this graph. 



1/T [sec- 1 ] 
800 




50 100 150 200 250 300 350 400 450 500 
No I a m [nm- 1 ] 

Fig. 2. The lifetime in trap of b = 2; The field is B = 0.1 G, 
B' = 6.1 kG/cm and ajjo = 82 nm. The lifetime is much 
shorter than that in Fig. 1. 



of F w (k) is roughly same as eq. (42). Hence, the increase 
of jl weakens the interaction between wfss and ns. The 
oscillation of the lifetime with a change of Nq appears in 
Fig. 2. 



4. Conclusion 

We have calculated the lifetime of condensate in a 
magnetic trap, where field's direction spatially changes 
and wfss, ns and sfss interact each other. Our investiga- 
tion has confirmed that the adiabatic approximation is 
applicable for quantum condensates. 

In the case that the condensate size is much smaller 
than the scale of change of field's direction (b ^> 1), the 
lifetime is considerably long and the condensate is stable. 
The sign of m does not concern the stability. When 
the population of trapped atoms changes, its distribution 
also changes and the lifetime oscillates as a function of 

No I ano ■ For certain distributions the condensates never 
decay with this process. 

The traps used in experiments are generally b 1 
and the decay of the condensate is not observable. The 
lifetime strongly depends on the energy of condensate 
and therefore on b. It is shortened as b decreases. If one 
uses a trap of b < 3, the decay of the condensate will 
be observed. Finally we give an example to observe this 
drastical change. When one sets the field of Bo = 1 G, 
B' = 6 kG/cm (b = 11) and compresses 50 atoms//xm 
of 87 Rb, the lifetime of this condensate is determind by 
three-body recombination as about 1 sec. After sweeping 
the bias field to 0.1 G adiabatically (b = 2), the lifetime 
becomes about 30 msec. This decay should be observed 
without disturbance by recombination of atoms. When 
the population of trapped atoms decreases, for instance 
10 atoms//im, the lifetime is shorten again to 3 msec. 
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